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Understanding of charge transport mechanisms in nanoscale structures is essential for the devel- 
opment of molecular electronic devices. Charge transport through ID molecular systems connected 
between two contacts is influenced by several parameters such as the electronic structure of the 
molecule and the presence of disorder and defects. In this work, we have modeled ID molecular 
wires connected between electrodes and systematically investigated the influence of both soliton for- 
mation and the presence of defects on properties such as the conductance and the density of states. 
Our numerical calculations have shown that the transport properties are highly sensitive to the po- 
sition of both solitons and defects. Interestingly, the introduction of a single defect in the molecular 
wire which divides it into two fragments both consisting of an odd number of sites creates a new 
conduction channel in the center of the band gap resulting in higher zero-bias conductance than 
for defect free systems. This phenomenon suggests alternative routes toward engineering molecular 
wires with enhanced conductance. 



I. INTRODUCTION 

The investigation of quantum transport through 
molecular systems has become an important research 
field in the last few decades X Progress in measuring and 
fabrication techniques 2 - has led to the continuous minia- 
turization of electronic devices, which have reached the 
point where quantum effects are important. For exam- 
ple, semiconductor devices have been reduced in size 
to the nanoscale 3 -^ and even to the atomic scale^^ In 
addition to the report of metal atomic wires* 7 - - — sta- 
ble and rigid carbon atomic chains have been reported 
recentl y 12 ' 13 A key idea behind the advances in the un- 
derstanding of charge transport through molecular sys- 
tems is based on the view proposed by Landauer, con- 
ductance is transmission 14 ^ 5 

Potential applications of molecular devices range from 
novel computer architecture o 16 ' 17 to chemical sensors^ 
and medical diagnostics^. Among the various molecu- 
lar devices, ID conductors such as molecular wires j22r— 
have been considered to be one of the most fundamental 
components for nanotechnology. Due to the reduction of 
size and dimensions of materials, ID systems show sensi- 
tive response to external field or intrisic characteristics, 
which can be exploited for the development of molecular 
devices such as biological sensors^ with high sensitivity. 

Polymers are one of the most promising materials 
for acting as ID conductors. Their applicability 2 ^ 
ranges from displays 2 ^— to thin film transistors^ 7 - 
photovoltaics^ and solar celrs i 29 ' 30 The compatibility of 
polymer materials with light-weight, mechanically flexi- 
ble plastic substrates and new fabrication methods make 
them possible candidates for future electronic devices. 

Intrinsic properties of molecules as well as external 
perturbations can strongly affect the transport proper- 
ties of the resulting low dimensional devices. For ex- 



ample, the conductance of ID polymers drastically de- 
pends on the concentration and position of impurities 
(dopant) 31 ' 32 or defects^ 3 - Other effects which are impor- 
tant include dimerization and the formation of solitonic 
defects. Hence, all of these properties should be consid- 
ered in the design of ID devices. 

Therefore, in order to facilitate the design of control- 
lable molecular devices, it is necessary to use a theoreti- 
cal approach that simplifies the monitoring of how these 
intrinsic features affect their transport properties. Such 
approach should allow one to gain understanding and 
deep insights into the physical origins of the behavior of 
these materials. The knowledge acquired through these 
models can be used for the interpretation and elucidation 
of experimental observations, as guidelines for the plan- 
ning of experiments, as well as in the design of molecular 
devices. 

In this paper we model molecular junctions where ID 
molecular wire is connected between two electrodes using 
a tight-binding approach. Then we systematically inves- 
tigate how the degree of dimerization, and the position 
of solitonic and binding defects affect the transport prop- 
erties of such ID systems. We calculate the electronic 
density of states (DOS) and the conductance using the 
Landauer model in terms of the equilibrium Green func- 
tions. 

This article is divided into the following sections. Sec- 
tion II briefly presents the theoretical framework and sec- 
tion III presents the effect of solitonic and binding de- 
fects on quantum transport through molecular systems, 
respectively. Finally section IV summarizes this article, 
stressing on the importance of the defect position on the 
linear conductor. 

The effect of an odd vs. an even number of sites in 
the DOS and in the conductance is summarized in Sup- 
plementary Information A, and the effect of dimerization 
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is summarized in Supplementary Information B. In ad- 
dition, the relationship between dimerization and length 
dependence of conductance is summarized in Supplemen- 
tary Information C. 



II. THEORETICAL FRAMEWORK 
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FIG. 1: Tight-binding description of the molecular junc- 
tion considered in this work. The on-site energy terms, 
a n ; n = 1, 2, • • • , N, L, R, are set to for the simplicity. The 
/3l/r are hopping integrals for left/right ID electrode. Vl/r 
and Vn.n+i; n = 1, 2, • ■ • , JV — 1 are transfer integrals for the 
molecule/electrode interface in the left /right side and hopping 
integral between nearest neighbors, respectively. The detail 
of theoretical framework is shown in Appendix A. 



We model a molecular junction by connecting a ID 
molecular wire between two electrodes. We describe 
this molecular system using a standard tight-binding 
Hamiltonian^— Hq that considers only tt orbitals. 
Here, the dimerization is represented by a sequence of 
alternating weak and strong bonds while the solitonic de- 
fect is described by a pair consecutive weak bonds. For 
clarity, we model the electrodes as ID wires with Hamil- 
tonians and Hr, which result in a fair representation 
of infinite reservoirs. Fig. [T] shows the molecular junc- 
tion considered in this work. The tunneling between the 
central molecule and the left and right electrodes is repre- 
sented by two matrix elements Vl and Vr. Note that the 
electrode parameters and tunneling amplitude reduced to 
into two independent parameters Tl and Tr,; the escape 
rates to the left and to the right electrode, respectively^ 
The details this of theoretical framework are described in 
Appendix A. 

In principle, we can expand this method to general 
molecular systems which include any type of orbitals. In 
this work, we simplified the model Hamiltonian by only 
considering 7r-orbitals since the contribution from other 
orbitals to the conductance near the Fermi energy is neg- 
ligible in the sp2 carbon systems which we addressed. 
However, in the case of the molecular wires including 
metal atoms ; 39 ' 40 where s-, p-, and/or d-orbitals play a 
role in charge transport, those orbitals have to be taken 
into account. 



III. RESULT AND DISCUSSION 

A. Influence of the position of soliton 

Dimerization and soliton formation are closely related 
processes which have been extensively studied. The 
dimerization has been examined the first in the context of 
metal-insulator transition where dimerization is known as 
Pcierls distortion^ while the soliton formation has been 
examined in the field of conducting ID polymers42 These 
effects strongly modify the conductance and other elec- 
tronic properties of ID polymers. Thus, such effects in 
single molecular wires connected between contacts also 
should be considered. In order to estimate the influ- 
ence of solitons on transport, we modeled a ID molec- 
ular system consisting of an odd number of sites cou- 
pled with electrodes, and examined the dependence of 
the transmission spectra, total DOS (TDOS) and local 
DOS (LDOS) on the position of a soliton. For simplicity 
we assumed that there is a single soliton on a molecu- 
lar wire. However, this approach can be also applied to 
larger numbers of solitons on molecular wires. 




FIG. 2: Schematic description of poly-acetylene-based (PA- 
based) molecular wires consisting of 19 sites, with single soli- 
tonic defect, connected between two ID electrodes. Depend- 
ing on the position of solitonic defect, 8 cases can be consid- 
ered. Transfer integrals for double bonds and single bonds 
are set to Vd = 1.2/3 and V s — 0.8/3, respectively. Two hop- 
ping parameters nearby the solitonic defect are set to 1.0/3. 
Coupling constants are set as Vl/r, = 0.8/3. Left/right ID 
electrodes are treated by Newns-Anderson mode l 43 ' 44 with 
«l/r = and /3 L/R = /3. 

Figure [2] shows modeled poly-acctylene-based (PA- 
bascd) molecular wires consisting of 19 sites connected 
between ID electrodes and allowing for possible 8 soli- 
ton positions along the molecular wires. Transfer in- 
tegrals associated with double bonds and single bonds 
are set to Vd = 1.2/3 and V s = 0.8/3, respectively. The 
two hopping parameters on either side of the soliton 



3 




FIG. 3: (a) Transmission spectra and (b) TDOS of the molecular system consisting of 19 sites connected between two leads 
parameterized by the positions of a soliton. 




FIG. 4: Surface plot of LDOS of molecular wire shown in 
Fig. [5] as a function of energy and site index. The localized 
peaks of LDOS at E — migrate along the wire in conjunction 
with the position of soliton in Fig. [2] 



site are set to 1.0/3. All on-site energies are set to 
a n = 0;n = 1,2, ■■• ,N. The coupling constants to 



the electrodes are set as Vl/r = 0.8/3. The ID elec- 
trodes are treated by the Newns-Anderson mode l 43 ' 44 
with «l/r = and /3l/r = Fermi energy is set to 
Ep = 0. 

Fig. [3] shows calculated transmission spectra and 
TDOS of the ID molecular systems with different soliton 
positions. Interestingly, depending on the position of the 
soliton, the transmission at the Fermi energy (E = E-p) 
changes drastically (See animation in Supplemental In- 
formation). When the soliton lies near an electrode, 
the resonances in the transmission spectra and TDOS 
at E = Ep broaden with lower peaks. In this case, in 
spite of having an odd number of sites, the features of 
the TDOS and transmission resemble those for the PA- 
based molecular junctions which have an even number of 
sites (See Fig. S4(a) or (c) in Supplemental Information). 
On the other hand, when a soliton lies in the middle of 
the molecular wire, the resonances in the transmission 
spectra and TDOS sharpen and get narrower. The value 
of the transmission at E = E-p approaches its theoretical 
maximum of 1.0. 

In order to elucidate the relationship between the po- 
sitions of a soliton and sharpness of resonant peaks, the 
distribution of LDOS along the molecular framework was 
calculated. Figure [4] shows the LDOS surface plot as a 
function of energy and site index. The LDOS surface plot 
has peaks near the Fermi energy (E = E-p) where the soli- 
ton is located. The LDOS peaks at E = Ep broaden and 
reduce in height when the soliton is located at the ter- 
minal of the molecular wire (Fig.[5Ja) and (h)), whereas 
the LDOS peaks at E = Ep sharpen and narrow when 
the soliton is located in the middle of molecule (Fig. [U[d) 
and (e)). 

This feature can be explained as follows. When a soli- 
ton lies at the terminal of a molecular wire, the local- 
ized state on the molecular wire strongly interacts with 
surface states of the lead broadening the LDOS peaks 
around E = Ep. This means an electron occupying the 
soliton state easily escapes back to the electrode where 
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the electron was injected from, so that the chance of the 
electron traveling through to the other electrode is re- 
duced resulting in low transmission at E = Ep. Mean- 
while, when a soliton lies at the center of a molecular 
wire, the localized state of the molecular wires interacts 
symmetrically and weakly with surface states of the two 
leads producing sharp LDOS peaks. 

This conclusion can also be understood in terms of the 
generalized symmetry condition4£~— The transmission 
peak close to the resonance is approximately described 
as T max (£ ~ Eo) = 4r L r R /(r L + I^) 2 , where E is the 
energy of the molecular orbital associated to the soliton. 
Tl/r follows a simple relation, r L/R oc |Vl/r||*l/r| 2 j 
where \&l/r is the orbital amplitude of the molecule at 
the left/right interface. Thus, symmetric coupling of the 
molecule with two electrodes, IY = Tr, which corre- 
sponds to the soliton localization at the center of the 
molecule in case 4 or 5 in Fig. [2j gives optimal trans- 
mission probabilities, while asymmetric coupling, as can 
be seen in case 1 or 8 in Fig. [2J gives lower transmission 
probabilities. 



B. Creation of transmission channels by a defect 

We have examined how the position of a soliton mod- 
ifies the conductance of molecular wires and seen that 
the conductance is greatly reduced when the soliton lies 
near electrodes. Here, we investigate the influence of the 
defect and demonstrate the counter-intuitive result that 
the introduction of a defect satisfying special conditions 
creates a transmission channel leading to higher conduc- 
tance than in defect free systems. 

In general, realistic molecular junctions will contain 
defects. The defects may induce disorder in the molecu- 
lar framework. In the presence of defects both the on-site 
energies and the hopping integrals are modified. There- 
fore, the transport properties of the molecular junctions 
will also be changed. Thus we need to consider the effect 
of defects on quantum transport. We will now examine 
the dependence of the transmission and DOS on the po- 
sition of the defect. Here we only consider fluctuations in 
the hopping integrals and introduce a defect as a reduced 
transfer integral in a molecular system with an even num- 
ber of sites. For simplicity we consider a single defect in 
the molecular system. 

Figure[5]shows the schematic description of ID molecu- 
lar wire consisting of 6 sites, with a single defect, coupled 
to two electrodes. Transfer integrals for nearest neighbors 
and coupling strength at the left and right interfaces are 
set to V n , n +i = P and Vl/r = /3/2, respectively The de- 
fect is introduced as a reduction of the transfer integral 
V between the two sites on either side of the defect to a 
value of (3/2. We considered 5 cases of different positions 
of a single defect on a molecular wire and investigated 
position-dependence of the transmission and DOS of the 
molecular wire. 

Figure [5] shows the transmission spectra and TDOS of 
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FIG. 5: Schematic description of a non-dimerized ID linear 
molecule consisting of 6 sites, with a single binding defect, 
connected between two electrodes. Depending on the posi- 
tions of the binding defect, 5 cases can be considered (Case 
1-5). Transfer integrals for nearest neighbors and coupling 
strength at left/right interface are set to V n , n +i = 1-0/3; n = 
1,2, ■•■ ,N — 1 and Vl/r = P/2, respectively. The transfer 
integral corresponding to the position of the binding defect is 
reduced to fi/2. 



the molecular wire shown in Fig. [SJ Depending on the 
position of the defect, the transport properties change 
dramatically particularly around E = E-p. Fig. [7] shows 
I-V curves of the molecular wire calculated in the linear 
response regime. Surprisingly, the transmission at E = 
Ep in cases 1, 3, and 5 is higher than for defect free 
systems (Fig. [5]). 

In order to further analyze this behaviour, we calcu- 
lated the LDOS as a function of energy and site index. 
Figure [8] shows the surface plot of LDOS as a function 
of electronic energy and site index. The position of the 
defect is illustrated by changing the colour of the LDOS 
curve from red to blue on either side. When a defect is 
introduced such that the molecular wire is divided into 
two fragments which have an odd number of sites, for ex- 
ample 3 + 3, the LDOS plot of each fragment, as shown 
in Fig. [5{c), resembles that of a linear chain having 3 
sites (N = 3) without defects, as shown in Fig. 3(a). 
The LDOS in Fig. (5Jc) shows non-negligible split peaks 
around E = Ep. Likewise when a defect is introduced 
such that the molecular wire is divided into fragments 
containing 1 + 5 or 5 + 1 sites, the LDOS plot of longer 
fragment resembles that of a linear chain having 5 sites 
(N = 5) without defects as can be seen in Fig. 3(b). 
The LDOS in Fig. [5{a) and (e) show non-negligible split 
peaks around E = Ep. This is why, despite including 
the single defect, these cases lead to higher transmission 
probabilities and TDOS around E = Ep than the defect 
free systems shown in Fig. [6] and Fig. [7] 

On the other hand, when the defect is introduced such 
that the molecular wire is divided into two fragments 
both having an even number of sites, for example 2 + 4, 
the LDOS plot of each fragment resembles that of a linear 
chain having an even number of sites without defects. 
The LDOS in Fig.^b) and (d) are small around E = Ep. 
Thus, these cases lead to low transmission probabilities 
and TDOS around E = Ep as shown in Fig. [5] 
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FIG. 6: (a) Transmission spectra and (b) TDOS of the molecular system consisting of 6 sites parameterized by the position of 
a defect. 




Voltage (V/p) 



FIG. 7: I-V curves of the molecular wires consisting of 6 sites 
with a single defect. The I-V curves of the molecular wires 
with a single defect for case 2 and 3 are shown in green and 
blue, respectively. As a reference, I-V curve for the defect free 
molecular wire is shown in black. 



IV. CONCLUSION 




In summary, in order to elucidate the factors control- 
ling charge transport through ID molecules connected 
between electrodes we have modeled molecular junctions 
using the tight-binding method. Then we calculated 
the the influence of the parameters including molecu- 
lar lengths, degree of dimcrization, odd-even effects, soli- 
ton formation, and defects on transport properties of the 
molecular junctions at equilibrium using the Landaucr 
approach combined with the Green's function formalism. 

The numerical calculations have shown that the trans- 
port properties at the Fermi energy (E-p = 0) dramati- 
cally change depending on the degree of dimerization and 
on whether the number of the sites is odd or even. It has 
been shown that dimcrization of the molecular wires as- 
sociated with Pcicrls distortion reduces the DOS at the 



FIG. 8: Surface plot of LDOS of molecular wires with a 
single defect as shown in Fig. [5] as a function of energy and 
site index. In each case LDOS plots are depicted in different 
colors across a defect. 



Fermi energy, leading to exponential decay of length de- 
pendence of conductance of the molecular wires. We 
also proved that the damping factor is closely related 
to the degree of dimerization of the molecular wires. In- 
terestingly, the longer chain system without dimcrization 
showed no decay in the conductance at the Fermi energy. 
In the case of molecular wires with an odd number of 
sites, the extra non-bonding state appears in the middle 
of the band gap contributing to high conductance at low 
bias. Additionally, in the molecular wires with an odd 
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number of sites, the conductance at the Fermi energy 
was independent of the length of the molecular wires. 

Our calculations also demonstrate that the transport 
properties are highly sensitive to the position of solitons 
and defects. When a soliton lies near one electrode it 
strongly interacts with surface states of that electrode 
giving rise to a low, broad peaks in the zero-bias con- 
ductance, while when a soliton lies in the middle of the 
molecular wire, the interaction of the localized state with 
the electrodes is weaker leading to high conductance at 
the Fermi energy with sharp resonant peak. Concerning 
on the calculation of the influence of defects, we obtained 
a counter-intuitive result that in some cases defects in- 
crease the zero-bias conductance. Our findings are that 
a defect which divides the molecular wire into two frag- 
ments both having an odd number of sites creates a new 
conduction channel and enhances the zero-bias conduc- 
tance. Since thermal fluctuation naturally creates this 
sort of defects that, at least for a while, strongly favors 
charge transfer, we might infer that this is one of the 
underlying mechanisms of transport at room tempera- 
ture. We believe that the presented systematic study 
could be predictive for a wider material {e.g. semicon- 
ductor/organic interfaces). Our results could guide the 
synthesis of novel molecular wires with enhanced conduc- 
tors. 
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Appendix A: Theoretical framework 

We model a molecular junction by connecting a ID 
molecular wire between two electrodes. For clarity we 
model the electrodes as ID wires. We describe this sys- 
tem using a tight-binding Hamiltonian, considering only 
7T orbitals. The tight-binding parameters used depend ex- 
plicitly on the interatomic distances in real space. Thus, 
we can introduce defects or disorder by directly con- 
trolling on-site energies (diagonal elements) and hopping 



integrals (off-diagonal elements) of the electronic tight- 
binding Hamiltonian matrices. 

Figure 1 shows the molecular junction considered in 
this work. Our starting point is to divide the entire ID 
system into three regions; left electrode, right electrode, 
and central molecular region. We write the electronic 
Hamiltonian as 

H = H L + V h +H c + V R + H H , 

where -J/l/Rj Hq, and Ul/r are respectively the Hamil- 
tonian matrices for the left/right electrodes, for the the 
central molecular region, and the matrix representing the 
coupling between the central molecule and the left /right 
electrodes. In this case, we considered only nearest neigh- 
bor interactions, so that the terms Vl/r are scalar val- 
ues. The matrix elements of Hq are given by [Hc]m,n = 
a n $m,n + V m>n {l - 5 min ), where a n and V m , n ar e on- 
site energies and hopping integral between m-th and n-th 
atomic orbitals, respectively. In this study, the electrodes 
are approximated as semi-infinite ID chains with a single- 
orbital per site, thus the Hamiltonians of the semi-infinite 
electrodes Hl/r have similar forms with Hq- For sim- 
plicity we only consider the interaction between nearest 
neighbors. The matrix elements of left/right electrodes 
are given by [.Hl/r] 

Hereafter we normalize all hopping integrals, /?l/r an d 
V m ^ n , in terms of (3 and set all on-site energies as 0, 
a n = 0;n= 1,2, ••• ,N,L,R. 

In the orthogonal basis the retarded Green function is 
given by 

G R (E) = [(E + \n)I -H c - E L (£) - Sr(S)]- 1 , 

where /, rj, and E L / R (-E) are the identity matrix, the 
positive infinitesimal, and the self energy in the left/right 
electrode, respectively. The self-energy term E l /r(-E) is 
defined as T, L/K (E) = V h/R g h/R {E)vl /R , where g h / R {E) 
is the surface Green function of left/right ID electrode. 
The left/right electrode consists of equally spaced ID 
sites and coupling between nearest neighbors is set as 
Pl/r = P- The surface Green function of the elec- 
trode consisting of equally spaced ID sites considering 
only nearest neighbor interaction can be analytically de- 
scribed by the Newns- Anderson mode l 43 ' 44 as g^/ R (E) = 
cxp (ifc)//?L/R whose derivation 4 ^ is summarized in Ap- 
pendix B. 

The conductance of a molecular junction in a low bias 
is estimated from Landauer's formula G = GqT(E), 
where Go = 2 (for spin)^- is the conductance quan- 
tum and T(E) is the electronic transmission probability. 
T(E) is obtained from the Fisher-Lee relation^: T(E) = 
Tr[G Ii (E)T L (E)G A {E)r R (E)}, where T h/K (E) is the 
broadening function defined as T L / R (E) = i[E L / R (£J) — 

Ej, R (E)]. The DOS is calculated from following equa- 
tion: 

DOS(£0 = Im [G(E) - G\E)} /2tt. 
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In this study we only consider coherent transport and 
ignore internal scattering effects such as inelastic trans- 
port or incoherent transport^ associated with electron- 
phonon coupling ! 52 ' 53 



Appendix B: Analytic solution of surface green 
function in semi-infinite ID electrode 



Here we derive the analytic solution of surface green 
function of a semi-infinite ID electrode consisting of 
equally spaced sites. Consider a semi-infinite ID elec- 
trode with lattice constant a as shown in Fig. [9] Each site 
has only one orbital whose position is given by x n = an 
and its on-site energy is a. We only consider the inter- 
action between nearest neighbors. The hopping integral 
and overlap between nearest neighbors are j3 and S, re- 
spectively Energy-dependent coupling between nearest 
neighbors is defined as f3 E = j3 — ES. 



s _ s 



a 



a t v2/ kJ %y * a 

P P fi (5 

n= 1 2 3 4 5 

M 



FIG. 9: Semi-infinite ID electrode consisting of equally spaced 
sites with lattice constant a. a is on-site energy. 5* and ft are 
hopping integral and overlap matrix between nearest neigh- 
bors, respectively. 

The wave function of semi-infinite electrode can be 
constructed from forward- and backward-propagating 
plane wave e ±lkx : 



1 



$(jfe) = -^-(c ikx -c~ ikx ). 
v2 



The bra and ket are defined as 

1 ( ,jkx g— ikx\\„\ 1 

1 /^—ikx rjkxWrr.X 1 



<Jr(fc)) = -L (e lkx - e- lkx )\x) = ^(2isinfcj;)|x) 
;*(fc)| = -^{e~ [kx - e ikx )\x) = 4=(-2isinfcc)(a:|. 



Thus, the wave function at site n is given by: 

(x n |*(fc)) = A=y2M^f{k)\x) = -^=(2ismkx r , 
V2 „ V2 



This function satisfies the boundary condition of semi- 
infinite ID electrode that the wave function on site n = 
is zero since ^(xq) — (xo\^(k)) = 0. 



The operator of Green's function is defined as 

|*(A;)X*(fc)| 



G R (E) = J2 



E + iS-e(k)' 



The free propagator of the electrode between sites n and 
n' is given by 



(G R (£))„,„, = £ 



{x n \*(k))(V(k)\x r . 
E + i5-e(k) 



Especially in the case of n = n! we get 



(G R (£))„,„ = £ 



2 sin kna 
E + iS- e(fe) ' 



We only need the matrix element at the terminal site 
since this is the only site where the molecular system 
couples, thus putting n = 1 gives 



(G R (£))u=E 



2 sin 2 fee 



E + i5-e(k) 



Hereafter we drop E from (G R (£'))i.i . Since the number 
of sites is very large thus the summation can be converted 
into integral -> 57 Jl„ /a dfc) leading to 

(g r )m = — r „ 2sin2fcfl „dfe. 



2tt J_, /a E + i5- s(k) 

Using E-k relation of ID electrode consisting of equally 
spaced sites (e(fe) = a + 2f3 E cosfea), we get 



(G% 



-7T j a 



sin ka 



7r J-K/a E + 16 — a — 2(3 E cosfea 



dfc. 



When we define 9 = ka, the variable of integral can be 
changed as 



(G R ) M 



sin 2 9 



1 



7T J_ n E + 15 - a - 2/3 E cos ( 



-d9. 



When we define z = e l9 , then we get (z — 1/z) 2 = 
—4 sin 2 9, (z + 1/z) = 2 008 0, and dz = izd9. Substi- 
tuting these expressions into (G^)i^i, we obtain 



(G R )i,i 



(z 2 - I) 2 



dz, 



4mp E J z 2 (z 2 -2pz + l + id) ' 



where p = -|^# = cos 9. The term z 2 — 2pz + 1 can be 
changed to (z — zi)(z — z-i). These Z\ and zi must satisfy 
following relations: Z\Z<x = \ and Z\ + Z2 = 2p. Note that 
Z\ = j- = e Wl satisfies these relations. Then (G^i^i is 
written as 



(G% 



(z 2 - I) 2 



4mfi E J z 2 ((z-zi)(z-z 2 ) + i<5) 



f(z)dz. 



dz 



4i7T^ i 



(Bl) 



We can solve this complex integral using residue theorem: 
(GR)l ' 1 = 4i^ 27ri ^ RCS[/(z)] ^- 



fc=0 



Figure [TU] shows the singular points of function f(z) 
on complex variable plane. The residual term for 
Rcs[/(z)] z=Z2 vanishes since the imaginary positive in- 
finitesimal term iS move the singular point at z 2 outside 
of the circle. Therefore, (G R )i,i is 



(G R )i,i 



2/3 J 



;(Re S [f(z)] z=Zo +Res[f(z)} z=Zl ), 



where 



Res[f(z)] z=Zl = [(z - zi)f(z)] z=Zl =zi- Z2, 




and 



Res[/(z)] 2 



(z 2 - 1) 



dz (z - zi)(z - z 2 ) 



zi + z 2 . 



Finally, we get the surface green function of ID electrode 
consisting of equally spaced sites as 



FIG. 10: Singularity points of the function f(z) on complex 
variable plane. The gray zone is surrounded by unit circle 
z = e lS . Because of the imaginary term 18 in Eq. (|B1[I . the 
singularity points, z\ and 22 move to inside and outside of the 
circle, respectively. 



where, V^ R = Vl/r, — -ESl/r is energy-dependent cou- 
pling strength between the electrode and molecule in 
non-orthogonal description, where Sl/r is the overlap be- 
tween the terminal site in the electrode and the molecule. 
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